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Coupled  energy-drift  and  force-balance  equations  that  contain  a  frictional  force  for  the  center-of-mass 
motion  of  electrons  are  derived  for  hot-electron  transport  under  a  strong  dc  electric  held.  The  frictional  force 
is  found  to  be  related  to  the  net  rate  of  phonon  emission,  which  takes  away  the  momentum  of  a  phonon  from 
an  electron  during  each  phonon-emission  event.  The  net  rate  of  phonon  emission  is  determined  by  the  Boltz¬ 
mann  scattering  equation,  which  depends  on  the  distribution  of  electrons  interacting  with  phonons.  The  work 
done  by  the  frictional  force  is  included  into  the  energy-drift  equation  for  the  electron-relative  scattering  motion 
and  is  found  to  increase  the  thermal  energy  of  the  electrons.  The  importance  of  the  hot-electron  effect  in  the 
energy-drift  term  under  a  strong  dc  held  is  demonstrated  in  reducing  the  field-dependent  drift  velocity  and 
mobility.  The  Doppler  shift  in  the  energy  conservation  of  scattering  electrons  interacting  with  impurities  and 
phonons  is  found  to  lead  to  an  anisotropic  distribution  of  electrons  in  the  momentum  space  along  the  held 
direction.  The  importance  of  this  anisotropic  distribution  is  demonstrated  through  a  comparison  with  the 
isotropic  energy-balance  equation,  from  which  we  hnd  that  dehning  a  state-independent  electron  temperature 
becomes  impossible.  To  the  leading  order,  the  energy-drift  equation  is  linearized  with  a  distribution  function  by 
expanding  it  into  a  Fokker-Planck-type  equation,  along  with  the  expansions  of  both  the  force-balance  equation 
and  the  Boltzmann  scattering  equation  for  hot  phonons. 
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I.  INTRODUCTION 

There  have  been  many  theories  proposed  over  the  years  to 
describe  the  transport  of  hot  electrons  through  a  solid  when 
the  electric  field  being  applied  to  the  system  is  large,  and  the 
current-voltage  characteristics  deviate  from  those  of  linear 
response  theory.  The  semiclassical  (regular)  Boltzmann 
transport  equation,  with  momentum-drift  included,  seems  to 
be  the  only  one  of  these  early  theories  amenable  to  practical 
use.  However,  even  the  regular  Boltzmann  transport  equation 
becomes  laborious  when  the  system  goes  beyond  the  linear- 
field  regime.  Frolich  and  Paranjape1  used  a  displaced  Max¬ 
wellian  distribution  to  describe  the  electron  transport  in  in¬ 
sulators  and  semiconductors  at  high  temperatures.  Later, 
Arai2  used  a  similar  model  to  describe  electron  transport  in 
metals  at  low  temperatures,  using  a  Fermi-Dirac  distribution 
instead  to  predict  a  finite  electron  temperature  under  an  ap¬ 
plied  electric  field,  even  when  the  lattice  temperature  went  to 
zero.  Soon  after  that.  Lei  and  Ting3  proposed  coupled  force- 
balance  and  energy-balance  equations  to  describe  electrons 
in  semiconductors  and  metals  at  both  low  and  high  tempera¬ 
tures  by  assuming  an  isotropic  quasi-thermal-equilibrium 
(Fermi-Dirac)  distribution  for  hot  electrons  with  a  tempera¬ 
ture  different  from  the  lattice  temperature.  Very  recently, 
Huang  et  al ,4  used  the  Boltzmann  scattering  equation  to  re¬ 
place  the  energy-balance  equation,  and  then  accurately  deter¬ 
mined  the  electron  temperature. 

It  is  well  known  that  the  drifting  of  electrons  under  a  dc 
field  can  be  treated  as  a  field-driven  center-of-mass  (collec¬ 
tive)  motion  of  many  electrons.5  The  scattering  of  electrons 
by  the  lattice  ions  or  impurities  within  the  lattice  can  then  be 
considered  as  a  relative  motion  within  the  center-of-mass 
frame.  A  spatially-uniform  external  field  will  only  couple  to 
the  center-of-mass  motion.  The  lattice  ions  and  the  ionized 
impurity  atoms  within  the  semiconductor  remain  stationary 


with  respect  to  the  moving  center  of  mass.  The  center-of- 
mass  motion  will  couple  to  the  relative  scattering  motion 
through  a  Doppler  shift  in  the  electron  energy  along  the  field 
direction.  This  can  be  modeled  as  a  frictional  force  acting 
between  the  drifting  electrons  and  the  stationary  lattice  and 
impurity  atoms.3  As  a  result,  the  motion  of  the  electrons 
under  a  dc  field  will  be  effectively  opposed  by  this  friction. 
Even  though  this  frictional  force  can  be  measured  classically, 
it  is  a  quantum-statistical  average  of  all  the  frictional  forces 
acting  on  all  the  electrons  in  the  relative  scattering  motion, 
and  thus  it  depends  on  the  distribution  of  electrons  in  differ¬ 
ent  quantum  states.  This  distribution  of  electrons  in  the  vari¬ 
ous  quantum  states  is  determined  by  the  Doppler-shift- 
modified  elastic  and  inelastic  scattering  of  the  electrons, 
including  phonon  and  impurity  scattering. 

The  most  straightforward  and  simplest  way  to  treat  a 
steady-state  distribution  of  electrons  is  by  introducing  a 
state-independent  electron  temperature  for  use  in  a  Fermi- 
Dirac  distribution  of  hot  electrons.  This  electron  temperature, 
different  from  the  lattice  temperature,  can  be  found  by  using 
an  isotropic  energy-balance  equation  for  the  relative  motion 
of  the  electrons.3  However,  this  simple  method  cannot  be 
used  when  an  external  field  is  present,  even  in  a  steady  state, 
since  in  that  case,  the  assumption  of  an  isotropic  Fermi-Dirac 
distribution  for  hot  electrons  cannot  be  justified.  Indeed,  the 
energy-balance  equation  itself  cannot  be  justified  when  the 
drift  velocity  is  large,  producing  a  Doppler  shift  that  is  com¬ 
parable  to  the  phonon  energy.  Finally,  even  though  the 
energy-balance  equation  can  include  the  screening  effect,  it 
excludes  the  effects  of  electron-electron  scattering  on  elec¬ 
tron  transport.  For  these  reasons,  we  introduce  a  more  rigor¬ 
ous  method  to  describe  the  relative  motion  of  hot  electrons 
undergoing  anisotropic  scattering  that  involves  the  use  of  the 
Boltzmann  scattering  equation,  with  the  addition  of  an 
energy-drift  term.  When  the  dc  electric  field  is  very  strong. 
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the  distribution  of  hot  electrons  is  far  from  equilibrium,  and 
the  definition  of  a  state-independent  electron  temperature  be¬ 
comes  impossible.  Therefore,  the  approximate  treatments3-6 
of  hot  phonons  are  no  longer  valid. 

The  motivation  of  the  current  paper  is  as  follows.  It  is 
known  that  the  regular  Boltzmann  transport  equation  in¬ 
cludes  carrier  drift  in  momentum  space,  but  does  not  apply  to 
hot-carrier  transport  under  a  strong  dc  field.  On  the  other 
hand,  the  energy-balance  equation  includes  hot-carrier  ef¬ 
fects,  but  the  isotropic  carrier  distribution  that  is  assumed 
does  not  include  any  anisotropic  momentum  dependence  in 
the  carrier  distribution  along  the  field  direction;  the  state- 
independent  electron  temperature  that  is  introduced  becomes 
unphysical  when  the  dc  field  is  very  strong,  and  the  condi¬ 
tions  are  far  from  equilibrium.  We  propose  coupled  energy- 
drift  and  force-balance  equations  that  can  be  applied  to  hot- 
carrier  transport  by  including  both  the  hot-carrier  effects  and 
the  anisotropic  momentum  dependence  in  the  carrier  distri¬ 
bution  along  the  field  direction,  without  introducing  a  state- 
independent  electron  temperature. 

The  energy-drift  and  force-balance  equations  proposed 
here  are  based  on  the  following  physical  considerations: 

(i)  Transport  of  carriers  under  a  dc  field  results  from  the 
center-of-mass  drifting  motion.  This  slow  motion  can  be 
treated  classically  after  a  quantum-statistical  average  has 
been  taken. 

(ii)  Interaction  between  moving  carriers  and  static  lattice 
ions  and  static  impurity  atoms  in  a  motion  relative  to  the 
center  of  mass  is  modified  by  a  Doppler  shift  in  the  energy  of 
the  moving  carriers  in  the  field  direction.  This  couples  the 
center-of-mass  motion  with  the  relative  motion. 

(iii)  The  relative  scattering  motion  of  the  carriers  becomes 
anisotropic  in  the  field  direction  due  to  the  Doppler  shift. 
This  contributes  to  a  nonzero  frictional  force  that  resists  the 
dc-driving  force.  This  frictional  force  contains  contributions 
from  impurities  and  phonons.  The  phonon-emission  process 
will  take  the  momentum  of  a  phonon  away  from  a  carrier, 
while  the  phonon-absorption  process  will  add  the  momentum 
of  a  phonon  to  a  carrier. 

(iv)  The  classical  center-of-mass  motion  can  be  described 
by  the  force-balance  equation  including  a  frictional  force, 
while  the  relative  motion  is  composed  of  many-particle 
quantum-scattering  events  that  can  be  described  by  the 
energy-drift  equation,  including  the  work  done  by  the  fric¬ 
tional  forces  to  increase  the  thermal  energy  of  the  carriers. 
The  increased  thermal  energy  of  hot  carriers  is  expected  to 
reduce  the  field-dependent  carrier  drift  velocity  and  mobility. 

In  this  paper,  we  have  generalized  our  theory  in  Ref.  4  by 
using  the  Boltzmann  scattering  equation  with  energy  drift 
included  and  by  using  momentum  dissipation  in  the  repre¬ 
sentation  of  the  phonon-induced  frictional  scattering  forces. 
This  allows  us  to  completely  eliminate  the  need  for  any  sort 
of  electron-  and  lattice-temperature  definitions,  and  allows  us 
to  describe  events  far  from  equilibrium.  This  approach  also 
brings  the  hot-electron  transport  formulation  closer  to  that  of 
the  semiconductor  Bloch  equations,7  which  should  allow  the 
coupling  of  these  two  formalisms  to  proceed  without  any 
problems  when  coherent  optical  interactions  are  included  in 
the  future. 

The  organization  of  the  paper  is  as  follows.  In  Sec.  II,  we 
introduce  our  model  and  theory,  and  derive  coupled  energy- 


drift  and  force-balance  equations  for  hot-carrier  transport  un¬ 
der  a  strong  dc  electric  field.  The  numerical  results  are  dis¬ 
played  in  Sec.  Ill  for  the  calculated  drift  velocities  and 
mobilities  as  functions  of  the  dc  field,  and  they  are  explained 
physically.  The  paper  is  concluded  in  Sec.  IV  along  with 
some  remarks. 

II.  THEORY 

The  motion  of  many  electrons  in  n-doped  bulk  semicon¬ 
ductors  can  be  separated  into  the  center-of-mass  and  relative 
scattering  motions.  The  center-of-mass  motion  of  electrons  is 
described  by  a  Newton-like  force-balance  equation  for  the 
drift  velocity  of  the  center  of  mass.  The  forces  in  this  equa¬ 
tion  contain  a  driving  force  from  a  dc  electric  field  and  a 
frictional  force  from  both  impurity  and  phonon  scattering. 
During  the  scattering  of  electrons  with  phonons,  each 
phonon-emission  process  takes  away  the  momentum  of  a 
phonon  from  an  electron.  Meanwhile,  each  phonon- 
absorption  process  adds  the  momentum  of  a  phonon  to  an 
electron.  The  relative  scattering  motion  of  electrons  should 
be  described  by  the  energy-drift  equation,  i.e.,  the  Boltzmann 
scattering  equation  without  momentum  drift  included,  but 
with  energy  drift  included.  The  energy  drift  is  due  to  the 
work  done  by  the  frictional  force  against  the  center-of-mass 
drift  motion  of  the  electrons.  This  work  increases  the  internal 
energy  (i.e.,  the  thermal  energy)  of  electrons,  causes  the  dis¬ 
tribution  of  electrons  to  deviate  from  an  equilibrium  one,  and 
reduces  the  mobility  of  hot  electrons  due  to  enhanced  scat¬ 
tering  with  phonons.  When  the  dc  field  is  weak,  the  phonons 
are  described  by  a  distribution  that  is  nearly  in  equilibrium. 
However,  when  the  dc  field  is  strong,  the  phonons  are  de¬ 
scribed  by  a  nonequilibrium  distribution.  For  this  latter  case, 
many  phonons  must  be  generated  to  balance  the  strong  dc- 
driving  force. 

In  order  to  make  this  paper  self-contained,  we  will  rewrite 
some  of  the  key  equations  from  Ref.  4  in  a  slightly  more 
general  form.  In  the  presence  of  a  spatially-uniform  dc  elec¬ 
tric  field  Edc,  the  Hamiltonian  of  many  interacting  electrons 
in  bulk  semiconductors  can  be  written  as 

K  =  ;p^2  Pi  +  2  - p — ~  <?2  n  ■  EAc 

2 m  i  i<j  4i76r0er|r,--ry| 

+  2  Uimp(r i  -  Ra)  -  2  ut  ■  V  ;.t/ion(ri  -  Rt),  (1) 

i,a  /,€ 

where  i=  1 ,2, . . .  ,Ne  is  the  index  of  Ne  electrons,  a 
=  1,2,...,  Na,  is  the  index  for  Na  impurity  atoms,  l 
=  1,2,...  ,N(  is  the  index  for  Nf-  lattice  ions,  r,  is  the  position 
vector  for  the  ;th  electron,  Ra  and  /p  are  the  position  vectors 
of  impurity  atoms  and  lattice  ions,  uf  represents  the  ion  dis¬ 
placement  from  the  thermal  equilibrium  position,  m  is  the 
effective  mass  of  electrons,  e0  is  the  dielectric  constant  in  the 
vacuum,  and  er  is  the  relative  dielectric  constant  of  host 
semiconductors.  The  single-electron  momentum  operator  is 

Pi=-ihV and  both  the  impurity  potential  U,mv(r  -  Ra)  and 
the  ion  potential  I/ion(r,-Rf)  are  included.  We  first  define  the 
center-of-mass  momentum  and  position  vectors  by 
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and  those  for  the  relative  motion  by 


Pi 


r\  =  r;-Rc. 


(2) 


(3) 


By  using  the  center-of-mass  and  relative  momentum  and  po¬ 
sition  vectors  defined  in  Eqs.  (2)  and  (3),  we  can  separate  the 
total  Hamiltonian,  including  the  Hamiltonians  of  electrons 
and  phonons,  into  one  center-of-mass  Hamiltonian  Hcm  and 
another  relative  Hamiltonian  7r!rel,  given  by 


H 


cm 


(pc)2  -  - 

—  j  -  NeeEdc  ■  R, 
2Nem 


(4) 


Ke\  =  2  Siflk,Aka+  2  h£lqJ)\J)qX 

k,a 

1  Y  el  -4  *t 

+  r,  ^  ^,ak+qo-ak’  -qtr'ak'  <T,ak<T 

^  r  ?t  t  €()€,-(]  V 
k,k  , a, a  C1  u  ' 1 

+  22  Cqk(bqx  + 

k,(T 


+  22  Ui(q)e“l  iRl'  Ra>dl+Saai (5) 

k,cr 


where  V  is  the  volume  of  the  system,  fiflqx  is  the  phonon 
energy  with  wave  number  q  for  mode  X  (totally  three 
modes),  sk=h2k2 12m*  is  the  kinetic  energy  of  electrons  with 
wave  number  k,  and  the  index  cr=  + 1  is  for  the  up-and  down- 
spin  states  of  electrons.  We  use  at  (dkn)  to  represent  the 

creation  (annihilation)  operator  of  electrons  and  h'--jX  (bqK)  to 
denote  the  creation  (annihilation)  operator  of  phonons.  U[{q) 


is  the  Fourier  transform  of  the  impurity  potential  U""p(r' 
+Rc-Ra),  and  CqK  is  the  electron-phonon  coupling  constant 
that  will  be  given  later  in  this  section.  The  coupling  of  the 
center-of-mass  and  relative  motions  [the  factor  exp(/c/-/C)] 
can  be  seen  from  the  impurity  and  phonon  parts  of  the  rela¬ 
tive  Hamiltonian  in  Eq.  (5). 

From  the  total  Hamiltonian,  we  derive  two  Heisenberg 
equations  for  the  center-of-mass  motion  of  electrons, 

=  I[pc,Kcm  +  nrc]\  =  NeeEdc  -  ;2  cqXqe^RC 

at  in 

q,K 

X(V  +  b-q^Pq  ~  '2  U^qe1^^ pq,  (6) 

q,a 

»=jRC=^[RC,Km  +  HrJ  =  ^,  (7) 

at  in  Nem 

where  pq=1<bcA\+qJ*klT  represents  the  density-wave  operator 
of  electrons.  Applying  a  quantum-statistical  average  ((■  ■  -))av 
to  Eqs.  (6)  and  (7)  and  defining  the  following  quantities: 


where  u0  =  ((Pc / Nem*))av  is  the  drift  velocity  related  to  the 
center-of-mass  momentum,  we  obtain  the  following  force- 
balance  equation: 

Npnjuo  =  NeeEAc  +  ^[hq]  +  Fp[w0] .  (9) 

Here  and  Tp[m0]  are  the  frictional  forces  resulting 

from  impurity  and  phonon  scattering,  respectively. 

For  the  relative  motion  of  electrons,  the  energy-drift  equa¬ 
tion  for  electrons  in  a  dc  field  can  be  derived  as8 


dn 
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k  )  [«] 


ds 


k)[n ] 


k,q  iyq^ 
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x  abs  I 

’U 


.(out) 


nb 


(10) 


where  nk  is  the  electron-distribution  function  in  the  k 
state,  and  the  energy-drift  term  -( |)[d/otejt]r„]{(Fi[wo] 

+  Lp[(70])  •  (7,)}.  due  to  the  work  done  by  the  frictional  force 
against  the  drift  motion,  has  been  included.  The  energy-drift 
term  physically  represents  the  rate  for  increasing  the  thermal 
energy  of  electrons  in  the  k  state.  The  work  done  by  the 
frictional  force  was  introduced  in  the  energy-balance 


equation3,8  to  determine  a  macroscopic  (state-independent) 
electron  temperature.  On  the  other  hand,  the  work  done  by 
the  frictional  force  is  introduced  here  to  microscopically  de¬ 
termine  the  thermal  effects  on  the  electron  distribution  under 
high  electric  fields,  without  having  to  define  an  electron  tem¬ 
perature,  which  is  unphysical  in  the  situation  described  here. 
The  introduction  of  this  additional  term  in  Eq.  (10)  is  crucial 
for  the  discussion  of  thermal  effects  on  hot-electron  trans- 
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port.  In  Eq.  (10),  {d/ dsk}ini  indicates  that  the  energy  deriva¬ 
tive  only  acts  upon  the  occupation  probability  nk  or  nk±^, 
hq  ■  u0  is  the  Doppler  shift  in  the  energy  of  moving  electrons 
along  the  field  direction,  h  a>qX  is  the  energy  of  acoustic 
phonons  in  the  X  mode  with  a  wave  vector  q,  hcoLO  is  the 
energy  of  the  longitudinal-optical  phonons,  U^q) 
=Zce2 /[e0er(q2 + 2,)V]  is  the  interaction  between  electrons 
and  impurities,  Zc  is  the  charge  number  of  ionized  donor 
atoms,  Q2=(e2/ €Ger){ni  1 Trh2){l>Tp'(T3n)U2  for  static  screen¬ 
ing,  and  cr3D=Ne/V  is  the  electron  concentration  in  the 
doped-host  semiconductors.  NqX(coqX)  is  the  nonequilibrium 
acoustic -phonon  distribution  which  satisfies  the  following 
Boltzmann  scattering  equation  for  hot  phonons: 

dNc-x(cooX)  , 

=  0^-xKJ  + 1]  - 

_  NqX{(oqX)  -  N0(hujqX/kBT) 

where  ©|™  and  ©|xs  are  the  rates  for  acoustic-phonon  emis¬ 
sion  and  absorption,  respectively,  t-x  is  the  relaxation  time 
for  acoustic  phonons  with  wave  vector  q  and  in  mode  X, 


N0U)  =  [exp(x)-l]  1  is  the  Bose-Einstein  function,  and  T  is 
the  initial  lattice  temperature.  For  optical  phonons,  the  non¬ 
equilibrium  distribution  Nq(coLO)  can  be  determined  in  a 
similar  way  by  the  following  equation: 

dN  A  coin)  , 

— ^  =  ©^-LO)  +  1]  -  ©f  ^ L0) 

Nq(ioLO)  -  N0(h(oLO/kBT) 

- ,  (12) 

where  0™  and  @^bs  are  the  rates  for  optical-phonon  emis¬ 
sion  and  absorption,  respectively,  and  t-  is  the  relaxation 
time  for  optical  phonons  with  wave  vector  q.  The  lack  of 
dependence  of  «L0  on  the  wave  number  q  will  alter  the 
momentum  exchange  between  the  electrons  and  the  optical 
phonons  during  a  scattering  process.  We  have  introduced  a 
relaxation-time  approximation  in  Eqs.  (11)  and  (12)  for  ad¬ 
ditional  phonon  scattering  other  than  the  electron-phonon 
scattering.  This  additional  phonon  scattering  includes  both 
phonon-phonon  interaction  and  boundary  scattering  of 
phonons. 

The  scattering-in  rate  for  electrons  in  the  final  k  state  is 


( '  )  2  77 .  ^  2  77 

Vvk=Na-r%  \Uiiq)\2[rik_q5(sk-  sk_q  +  hq  ■  u0)  +  nZ+qS(sk- sk+q-hq  ■  u0)]+  \CqX\2{nk_^N^x((oqX)S{ek- sk_q  -  hcoqX 

11  i  n  ?,x 

2  77 *  - 

+  hq  ■  u0)  +  nk+^NqX(wqX)  +  1]  S(ek  -  ek+q  +  huqX-  hq  •  w0)}  +  \Cq\2{nk_qNq{uBO)  S(sk  -  ek_q  -  hu)L0  +  hq  ■  u0) 

n  q 


277i 


+  nk+^.Nq((oLO)  +  l]^(eJfc  -  ek+q  +  h(oLO  -  hq  •  w0)}  +  —  2 


„2  \  2 


h  .  \  e0e,q2V 


(1  -  nk,)nk_qnk,+q8(ek  +  sk,  -  sk_q  -  ek,+q),  (13) 


where  the  last  term  in  Eq.  (13)  is  due  to  Coulomb  scattering,  and  the  scattering-out  rate  for  electrons  in  the  initial  k  state  is 


Wiat)=Na^\uiqm 

h  - 


277i 


■  nk+q)  S(Ek+q  -ek  +  hq  •  Mq)  +  (1  -  nk_q)  S(sk_q  -sk-  hq  ■  u0)]  +  —  2  |  C?x|2{(  1  -  nk+q) 


q,\ 


2  77  - 

XNqX((oqX)S{ek+q  -  sk  -  haqX  +  hq  ■  u0)  +  (1  -  nk_q)[NqX(<oqX)  +  1  \S(sk_q  -ek  +  hcoqX  -  hq  •  w0)}  +  —  2  I Cq 

n  q 


x{(l  -nk+q)Nq(ML0)S(ek+q-  Ek-h(oLO+hq  ■  u0)  +  (1  -  nk_q)[Nq((oLO)  +  1  ] S(sk_q  -  sk  +  h(ol  0  -  hq  ■  il0)\ 
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2 

nk,{\  -  nk_q)(\  -  nk,+q)S(sk_q+  ek,+q 


~  ek~  eA :')■ 


(14) 


For  the  scattering-in  rate,  the  first  term  in  Eq.  (13)  is  dia- 
grammatically  represented  by  the  top  panel  in  Fig.  1.  The 
second  and  third  terms  in  Eq.  (13)  are  diagrammatic  ally  rep¬ 


resented  by  the  top  panel  in  Fig.  2.  On  the  other  hand,  the 
first  term  in  Eq.  (14)  for  the  scattering-out  rate  is  diagram- 
matically  represented  by  the  bottom  panel  in  Fig.  1.  The 
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FIG.  1.  Diagrams  for  scattering-in  (top  panel)  and  scattering-out 
(bottom  panel)  rates  of  moving  electrons  with  static-impurity  at¬ 
oms.  The  solid  lines  with  arrows  represent  the  electron  states  with 
wave  vectors  indicated.  The  dashed  lines  with  crosses  represent  the 
interactions  of  electrons  with  impurities  with  the  wave  vector  q. 
The  circled  wave  vector  k  of  electrons  denotes  the  initial  state  for 
the  scattering-out  rate  and  the  final  state  for  the  scattering-in  rate.  u0 
is  the  electron  drift  velocity  with  its  direction  indicated  by  a  hollow 
arrow  in  the  dc-field  direction.  The  conditions  for  energy  conserva¬ 
tion  in  each  case  are  also  shown  among  the  electron  kinetic  energies 
sk  and  ek±q,  and  the  Doppler  shift  hq-uQ. 

second  and  third  terms  in  Eq.  (14)  are  diagrammatic  ally  rep¬ 
resented  by  the  bottom  panel  in  Fig.  2. 

The  emission  rate  for  acoustic  phonons  in  the  X.  mode 
with  a  wave  vector  q  due  to  the  interaction  of  phonons  with 
electrons  is 


4tt 


=  E  =  —\ Cj2E  -  "*) 


J  k,q\  ^ 


4  TT 

kt-iU  -  »£) 

k  k 

X  S(ek  -  sk_q  -hcoqX  +  hq  ■  u0) 

+  nk{  1  -  nk+q)  S(ek+q  -sk-  h(oqX  +  hq  ■  u0)}. 


For  optical  phonons,  we  have 
477, 


(16) 


0™  =  E  lc9|2E  {nk+q(  1  -  nk) 

k  n  k 

X  S(ek  -  sk+q  +  hwL0  -  hq  ■  u0) 

+  nk(  1  -  nk_q)  S{sk_q  -ek  +  hwL0  -  hq  ■  u0)}, 


and 


(17) 


477, 


®fs  =  E  ®f-j  =  -J-\ C/E  {nk-q(  1  -  nk) 

k  ’  n  k 

x  -  ek -q  -  h(x)LO  +  hq  ■  il0) 

+  nk(]  -  nk+q )  -ek-  h(oLO  +  hq  ■  u0)}. 

(18) 

For  the  phonon-emission  rate  in  Eq.  (15),  its  diagrammatical 
representation  is  shown  by  the  top  panel  in  Fig.  3.  On  the 
other  hand,  the  diagrammatical  representation  of  the  phonon- 
absorption  rate  in  Eq.  (16)  is  displayed  by  the  bottom  panel 
in  Fig.  3.  In  Eqs.  (13)— (18),  we  have  defined9 


|CJ2  = 


h 


2ptVokq(  _ 


D~q~+—(ehu) 


<r 


q2  +  Q] 


(19) 


x  -  ek+q  +  hwqX-hq-  u0 ) 

+  nk(  1  -  nk-q )  S(Bk-q  -Zk+hMq\-  ^  '  «o)}^ 


(15) 


and  the  absorption  rate  for  acoustic  phonons  in  the  A.  mode 
with  a  wave  vector  q  is 


FIG.  2.  Diagrams  for  scattering-in  (top  panel)  and  scattering-out 
(bottom  panel)  rates  of  moving  electrons  with  static-lattice  ions. 
The  solid  lines  with  arrows  represent  the  electron  states  with  wave 
vectors  indicated.  The  dashed  lines  with  arrows  represent  the  inter¬ 
actions  of  carriers  with  phonons  with  the  wave  vector  q.  The  circled 
wave  vector  k  of  electrons  denotes  the  initial  state  for  the  scattering- 
out  rate  and  the  final  state  for  the  scattering-in  rate.  klqX=(oqX  or 
<uL0  denotes  the  phonon  frequency  and  q-ug  is  the  Doppler  shift. 


IQ, |2 


h  13 
2piVw?(  64 


2 


(20) 


for  the  couplings  between  the  electrons  and  the  longitudinal 
(\  =  L)  and  transverse  (X  =  ?)  acoustic  phonons  with  wqX 
=  qsx  in  the  Debye  model,  where  sx  is  the  sound  velocity  of 
acoustic  phonons  in  the  \  mode,  p,  is  the  ion-mass  density,  D 
is  the  deformation  potential,  and  /z  14  is  the  piezoelectric  con¬ 
stant.  For  optical  phonons,  we  have1 


FIG.  3.  Diagrams  for  phonon-emission  (top  panel)  and  phonon- 
absorption  (bottom  panel)  rates  due  to  interaction  of  moving  elec¬ 
trons  with  static-lattice  ions.  The  solid  lines  with  arrows  represent 
the  electron  states  with  wave  vectors  indicated.  The  dashed  lines 
with  arrows  represent  the  interactions  of  carriers  with  phonons  with 
the  wave  vector  q.  flqX  =  toqX  or  io^q  denotes  the  phonon  frequency 
and  q-u0  is  the  Doppler  shift. 
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-»  -> 


FIG.  4.  Illustration  for  orientations  of  electron  ( k )  and  phonon 
(i q)  wave  vectors.  The  applied  electric  field  Vc  and  electron  drift 
velocity  ii0  lie  in  the  ( x )  direction.  The  electron  wave  vector  k  can 
be  decomposed  as  a  parallel  component  (kx)  (in  the  x  direction)  and 
a  perpendicular  component  k±  (in  the  y  direction  due  to  rotational 
symmetry).  The  phonon  wave  vector  q  can  also  be  decomposed  in 
the  same  way.  The  angle  between  q ,  and  the  y  direction  within  the 
y-z  plane  is  denoted  by  <fi ,  and  the  angle  between  k  and  the  x 
direction  within  the  x-y  plane  is  represented  by  0. 


.  n  flG)\  n  /  1  1  \  €2 


(21) 


where  es  and  e,x  are  the  static  and  high-frequency  dielectric 
constants  of  host  semiconductors. 

By  assuming  a  dc  electric  field  along  the  x  direction,  as 
shown  in  Fig.  4,  the  force-balance  equation  for  the  center-of- 
mass  motion  of  transported  electrons  is  cast  into  the  form  of 


Nem*^f=NeeEic  +  Fx[u0\, 
at 


(22) 


where  the  quantum  statistically  averaged  frictional  force  is 
found  to  be  Fx[uq]  =  F'^uq\  +  F^"[uq\  with3-8 

F'x[u0]  =  ~  y-2  Hx(nk+q-  n^NilUiiq)]2 

k,q 

X  8(ek+q  -  ek  +  hqxUo) ,  (23) 


Fx\uo\  =  -  2  +  1]  -  ®tjjxNfa((QqX)} 

k,q,X 

-  2  hqJ^f^N^  o)  +  1]  -  0^-(wLO)}. 


F'x[u0]  and  are  due  to  impurity  and  phonon  scatter¬ 

ing,  respectively.  Because  of  the  Doppler  shift  in  LO-phonon 
energy,  i.e.,  (oLO-fiq- u0  for  phonon  absorption  and  emission 
rates  in  Eqs.  (17)  and  (18),  the  LO-phonon  contribution  to 
the  frictional  force  in  Eq.  (24)  becomes  nonzero.  In  addition, 
Eq.  (24)  reduces  to  -Nem*u0l  rph  in  the  leading  order  of 
small  u0  under  a  weak  electric  field,  where  rph  can  be  viewed 
as  a  momentum-relaxation  time  of  electrons  from  phonon 
scattering.  The  diagrammatical  representation  of  Fvx\uq\  in 
Eq.  (24)  can  be  seen  from  Fig.  5  for  the  phonon-emission 
event  (left  panel)  and  the  phonon-absorption  event  (right 


[NqUlJ 

+{Dqrq-Uo] 


[iV  q-'T] 


FIG.  5.  Diagrams  for  changes  of  the  electron  momentum  due  to 
phonon  emission  (left  panel)  and  phonon  absorption  (right  panel). 
For  each  phonon-emission  event  with  a  rate  and  an  occupation 
factor  Nqx+ 1,  the  electron  wave  number  is  reduced  from  the  initial 
value  kx  to  the  final  value  kx-qx  in  the  field  direction  by  the  phonon 
wave  number  q„  as  indicated  by  a  hollow  arrow  in  the  left  panel. 
For  each  phonon-absorption  event  with  a  rate  0“^  and  an  occupa¬ 
tion  factor  V?x,  the  electron  wave  number  is  increased  from  the 
initial  value  kx  to  the  final  value  kx+qx  in  the  field  direction  by  the 
phonon  wave  number  qx,  as  indicated  by  a  hollow  arrow  in  the  right 
panel.  The  total  change  of  the  electron  wave  vector  in  the  direction 
perpendicular  to  the  field  is  zero,  since  the  individual  contributions 
cancel  each  other.  denotes  the  phonon  frequency  and  q-ug  is 
the  Doppler  shift. 


panel).  The  diagrammatical  representation  of  F‘r[n0]  in  Eq. 
(23)  can  be  found  from  Fig.  6. 

When  ek,  s k±q> hq xu0,  ti ojqk  is  assumed  for  high  densities 
of  electrons,  we  can  expand  the  energy-drift  equation  up  to 
the  second  order  by  introducing  a  continuous  distribution 
function  f(sk±,skJ  =  p(e*)n*=[(2m*)3/2(eA)1/2/(2TT2^.3)]nt, 
with  ek=Bk  +ek i  Bk  =h2k2L/2m*,  ek  =h2k2x/2m* ,  k 
=  (kx,k1),  and  q  =  (qx,q1).  When  the  acoustic -phonon  fre¬ 
quency  (dominant  phonon  modes  at  low  temperatures)  and 
the  Doppler  shift  are  smaller  than  the  plasma  frequency  of 
electrons,  the  static-screening  model  can  be  justified.  For 
steady-state  electron  transport,  the  Coulomb  scattering  effect 
can  be  approximated  by  a  homogeneous  level  broadening 
related  to  the  lifetimes  of  quasiparticles.  By  using  the  above 
assumptions,  the  energy-drift  equation  leads  to  the  following 
linearized  Fokker-Planck-type  equation,  with  respect  to  the 
distribution  function, 


nk(1-nk.q)-nk-q(1-nk) 


-> 


[u0] 
iT  tq 

FIG.  6.  Diagrams  for  changes  of  the  electron  momentum  due  to 
scattering  with  impurities.  For  each  scattering  event,  the  electron 
wave  number  is  reduced  from  the  initial  value  kx  to  the  final  value 
kx-qx  in  the  field  direction  by  the  wave  number  qx  of  an  impurity 
atom,  as  indicated  by  a  hollow  arrow  at  the  bottom  of  the  figure.  u0 
is  the  electron  drift  velocity  with  its  direction  indicated  by  a  hollow 
arrow,  and  _/Va|l/i(g)|2  represents  the  strength  of  impurity  scattering. 
nk(l-nk-j)-nk-j(l-nk)  represents  the  phase-space  filling  effect 
(Pauli  exclusion)  when  the  electron  state  k  is  the  initial  or  final  state 
in  two  opposite  scattering  events.  The  total  change  of  the  electron 
wave  vector  in  the  direction  perpendicular  to  the  field  is  zero,  since 
the  individual  contributions  cancel  each  other. 
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jf(Bk±’Ek)  =at (Bk±,ek)f(ek±,ek)  -  [Vt(s*±,£*a) 
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which  is  subjected  to  the  conservation  of  the  total  number  of 
electrons. 
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r+oo 

dkLk_L 

Jo  J-c 


dk. 


f(Bk±,ek) 


p(ek ,  +  e* ) 


(26) 


The  spontaneous  phonon-emission  rate  is 


x  d(ek  ~  sk-q  ~  AcoqX  +  Aqxii0)  +  [NqX(a>qX)  +  1] 
X  S(ek  -  ek+q  +  h(oqX  -  hqxu0)} 
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(30) 
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The  thermal  (T)-and  dc-field  (F)-induced  energy  transfer 
rates  are 

Vj(BkrBk)  =  ~  27rA^iM02  qx\UXq)\2[S(ek -  ek_q  +  fiqxu0) 

q 

-  -  ek+q  -  HxU o)]  +  2tt2  (<0qk  -  qxu0 ) 

q,\ 

x  I C9>|2{^(<u?x)  S(ek  -  sk_q  -h(»qX  +  hqxu0) 
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The  thermal  (T)-and  dc-field  (F)-induced  energy-diffusion 
rates  are 
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Similar  to  the  second-order  expansion  of  the  energy-drift 
equation,  we  can  expand  the  frictional  force  up  to  the  first 
order. 


Fx\uq\  =  277^2  ~r^Ni(l2xuo\Uiiq)\ls^sk+q  ~ek  +  Aqxu0) 


P(e*) 


X 


k,q 

d  d 

+ 

dek  dsk 


f(sk .  >  ek) 


-  4tt2  ?.v|C9,\|2E - -[<5(e*  -  Bk+q 

q,X  k  P(£k) 

+  h(jjqk  -  hqxu0)  +  S(ek_q  -ek  +  ho)qX  -  hqxu0)] 

+  477*2  ^|C?,x|22  — -t{(9.v«o-  ^vJIXa^x) 

«,X  k 

+  lMs*  -  ek+q  +  ha)qX  -  hqxu0) 

+  ( qxuo  +  w9x)^x(w<yx) 

x  s(ek  -  ek+  -  hoj  k  -  hqxu0)}  - - 

dskl 


dek 


f(,Bk  ,Bk )  -  4rr2  ?.v|c9|22  — - 
L  X  a  i  p(e*) 


X  [<5(eA.  -  ek+q  +  Aml0  -  hqxu0)  +  S(ek_q  -  ek  +  fi&)L0 
-  Aqxu0)] 


+  4tt *2  ?JC,|22  “7~: 

9  k  P 


{( qXu0  ~  wLo)[-^o(wLo) 


+  l]^(e*  -  Bk+q  +  h(oLO  -  hqxu0 ) 

+  (qxu0  +  (olo)Nj((olo)  S(ek  -  ek+q  -  h(oLO  -  hqxu0)} 


X 


d  d 

—  "I"  2 

dek  dek 


f(Bk±,Bk). 


(32) 
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By  using  the  same  expansion  technique,  the  emission  and 
absorption  rates  of  acoustic  phonons  can  be  rewritten  as 

4tj-  f(.Bk  ’^k') 

©'“  =  —  |C^|22  _,X_  /  [<5 {ek  -  ek+q  +  h(0qX  -  hqxu0) 


h 


p(e*) 


+  8{ek-q  -  sk  +  hokqk  -  hqxu0 )]  +  4tt|C^|2(w9X  -  qxu0) 

X  Z  "TT  <5{ei  -  Sk+q  +  hvqk  ~  hqxu 0) 
k  Pk) 


X 


P  d 
—  +  — 
ds ds k 


f(skL’Sk)’ 


(33) 


^abs  4^lr  pv-^K’^L  r  r  , 

©9-x  =  -4  ^(e  ^  L^k-s^  -  ficoqX  +  hqxu0) 

+  8{ek+q  -sk-  fuoqX  +  hqxuQ)]  -  4tt|C9X|2(w9X  -  qxu0) 


1 

xZ~ y 

k  P(Bk) 


<5k  -  ek  -  h<o  k  +  hqxuQ) 


X 


P  d 

+ 

dek  dek 


f(Bk±,ek)- 


(34) 


For  optical  phonons,  we  have  similar  results. 


jaem.i ”\r  l2Y 

«  h  '  J  f  P(Bk) 


[<5k  -  ek+q  +  huLO  -  hqxu0) 


+  S[sk_q  -  ek  +  fiok10  -  hqxu0 )]  +  4tt|C9|2(wl0  -  qxu0) 

X  2  TT  <5k  -  Sk+q  +  t>  «LO  -  Hxuo) 
k  Pk) 


X 


P  d 
~  +  — 

dek  dek 


f(Bk±,ek), 


(35) 


©f s  =  TrlC/Z  y  /'  [<5k - 

11  k  P^ek' 

+  %Bk+q  -ek-  hw LO  +  hqxu 0)]  -  47r|C?p(wLO  -  qxu0) 
X  2  T  \  S(ek  ~  ek -q  ~  +  ftPxuo) 


Pk) 


X 


P  P 

“  Z 

dek  dek 


f(Bk±,ek). 


(36) 


By  defining  the  following  three  dimensionless  functions  after 
the  angle  integration  over  cf). 


rlir 

M)  = 

J  0 


dcf)  I 

27t| 


Jq/tt 


75  +  kv  ±  2  kxqx  +  q]_±2k±q1_  cos  4>  ±  2  qxu0f 


em 


.{?})=  f 

J  c 


_  r^l. 


yo/w 


i/kW,  {<?})=  I  T“ 


0  277  [  To  +  k*  ±  +  P1  ±  2£rtf±  cos  (f>±2(qxu0  -  qsx)f  J 

yo/rr 


■■ri 


0  277  I  75  +  k*  ±  2*a  +  <7±  ±  2k±  <3T  cos  4>  ±  2(qxu0  -  O)]2 


(37) 


(38) 


(39) 


where  {£}=(£.,., k±),  k±  =  k±/kp  kx=kx/kj,  k=(PL+k~)m,  {q} 
=  (qx,qA),  q±  =  q±lkf,  qx=qjkf,  q=(q2L+q2c)112,  u0 
=m*u0/hkf,  sx=m*sx/hkp  fi=m*flL0/(M2),  7o=7o ^E/> 

=  (3  t^ctid)  1/3>  Bj  =  h2kj/2m*,  and  y0  is  homogeneous  level 
broadening  due  to  Coulomb  scattering,  the  Fokker-Planck- 
type  equation  can  be  cast  into  a  dimensionless  form, 

d-_  _ 

-^f(ek±,Bk)  =  AT(ek±,ek)f(ekl,ek)  -  [VT(ey,e^) 

+  Vf(Bk]  ’kx)] 

+  [DT(Ek^,sk)  +  D¥{sk  ]  ,ep] 


d  d 

z  z 

ds ds k 


f(Ek±,ek) 


X 


P2  J?_ 
del  del 


f(ek±’Bk)> 


(40) 


which  is  subjected  to  the  constraint  [see  Eq.  (26)], 


r+  oo  r+oo 

dkj± J  dkx(p1+pxr,/2f(ski,sk)=—2,  (4D 

where  T=eft/h,  f(ek±,ek)  =  {ef/3Tra3D)f{ek±,ek)  =  [(Ek± 
+  eky/2/2n2\nk,  ek  =ek  / Ep  ek  =ek  / £y.  The  expansion  co¬ 
efficients  in  Eq.  (40)  can  be  found  from  Eqs.  (A1)-(A5)  in 
Appendix.  In  a  similar  fashion,  we  can  also  express  the  dy¬ 
namical  equation  for  hot  phonons  into  a  dimensionless  form 
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dNgX(2qsk) 

dt 


=  ®°£[N?x(2qsk)  +  1]  -  ®fkNjk(2qsk) 


N^(2qsx)-N0(2qsk/T) 


(42) 


and  (43),  the  expansions  for  the  dimensionless  emission  and 
absorption  rates  of  phonons  are  given  by  Eqs.  (A9)-(A12)  in 
the  Appendix.  In  addition,  the  force-balance  equation  can 
also  be  cast  into  a  dimensionless  form. 


dNp.VL) 

dt 


=  0|m[A^(2n)  +  1]  -  ©f Nqi2il) 


du  o 
dt 


-  £dc  + -F.v[moL 


(44) 


Nfr(m)-N0(2Cl/D 


(43) 


where  T=kBT/ej,  Tq\=EfTqk/h  and  7--=£jT-//L  In  Eqs.  (42) 


where  Edc=eEdc/ kfEp  and  the  dimensionless  frictional  force 
takes  the  form  of 


-  r  n  Fx[u0]  n-do  f  -  -  f 
FJw0]  =  =  — —  dk±k l 

NeSfkf  2^  Jo  J- 

d 


r+  00  r+ 

dq^q± 

J  0  J-c 


dkLkL  dkx(F±  +  Fx)  1/2  dqLq±  dqxq]\U{{q)\2  F{{k},{q}) 


dEu 


r+ oo  r+co  r+co  r+co 

f(B/c±,%)  -  dkLkL  dkx(k \  +  k^)-mf(ek±,ek)^j  dq±q±  dqxqx\C^\\^\{k},{q}) 

^  J  o  J  —oo  x.  *  o  J  — 

r+co  r+oo  r+co  -poo 

4_)({^}>®)]  +  3rr  dk  LkL  dkx(F±  +  ^)_1/22  dq±q±  dqxqx{qxu0  -  ^sJlC^J2  X  {[A^x(2pTJ 

J  0  —oo  X  J  0  —oo 


l]£x+)(U},{?})  +  Nfr(2qsk)i{  }({k}, {<?})} 


d  d 


377  r+co  -  -  r 
—  dk±k± 

^  J  0  j  -  c 


f(ski,sk)  -  —  J  dk±k±  j  dkjk 


12 

_L 


+  %rmf{ek 


r+co  r-t 

,  £k)  dq  i  q 
J  ()  J 


dek ^  de 

r+co  r- + 

dqLqL  |  dqxqx\Cg\2[i//+\{k},{q})  +  if^~\{k},{q})]  +  37rJ  dk±kLJ  dkx(k 


12 

_L 


r+co  r- 1 

+  £?r1/2  dqxqx(qxu0  -  n)|Q|2{[AU(20)  +  1  ]^+)({k},{q})  +  ^(2^)  (/^(W,  {(/})} 

J  o  J  -< 


X 


d  d 
Z  T-  _ 
dek 


/(«*,- e^’ 


(45) 


where  n=N.JV,  |C^J2,  |C^|2  can  be  found  from  Eqs. 
(A6)-(A8)  for  \=t,  and  |t7i(^)|2  can  be  found  right  below 
Eq.  (A5). 

III.  NUMERICAL  RESULTS 

For  our  numerical  calculations,  we  have  chosen  GaAs  as 
an  example  for  the  host  semiconductor.  For  GaAs,  we  have 
taken  the  following  parameters:  (r3D=  1  X  1018  cm-3,  n  =  1 
X  1015  cm-3,  m*/m0= 0.067,  Zc=  1,  ev=  1 3,  e„=  1 1 ,  e,.=  12, 
s^=5.14X  10s  cm/s,  s(=3.04X  105  cm/s,  D  =  - 9.3  eV,  hl4 
=  1.2X  107  V/cm,  p,  =  5.3  g/cm3,  y0=5  meV  and  T-k=T-t 
=  Tq= 3.5  ps.  The  other  parameters,  lattice  temperature  T  and 
electric  field  Edc,  will  be  given  in  the  figure  captions.  The 
electric  field  is  assumed  to  be  in  the  x  direction. 


Figure  7  displays  the  time  evolution  of  the  calculated  drift 
velocity  u0  as  it  gradually  reaches  a  steady  state.  This  time- 
evolution  method  has  been  used  previously8  to  seek  a  stable 
steady-state  solution.  In  this  case,  the  dc  field  is  turned  on 
right  after  t= 0.  From  Fig.  7,  the  steady  state  is  reached  with 
Uq  on  the  order  of  107  cm/s  after  t>  1.2  ps,  when  T=30  K 
and  Fdc=0.75  kV/cm. 

In  order  to  elucidate  the  physics,  we  show  the  calculated 
distribution  functions  of  electrons  (ef/kj)f(ek^,ek)  in  Fig.  8 
at  T=30  K  and  Edc=0.75  kV/cm  as  functions  of  ek  / Sj  with 
k±  =  0  in  Fig.  8(a)  and  skJsf  with  kx=  0  in  Fig.  8(b).  For  the 
purpose  of  comparison,  the  equilibrium  Fermi-Dirac  distri¬ 
bution  function  is  also  shown  in  Fig.  8  by  the  dashed-dotted 
curve  with  the  symbol  □.  From  Fig.  8(a)  it  is  easy  to  see  that 
electrons  are  moved  from  kx> 0  states  (dashed  curve  with 
the  symbol  A)  to  kx< 0  states  (solid  curve  with  the  symbol 
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FIG.  7.  Time  evolution  of  calculated  drift  velocity  u0  as  a  func¬ 
tion  of  time  t  at  T=  30  K  and  Eic=0.15  kV/cm  for  hot  electrons. 
The  dc  field  is  turned  on  right  after  t= 0.  u0  linearly  increases  with 
t  initially,  and  a  steady  state  is  finally  reached  after  t>  1.2  ps. 

O)  in  the  field  direction,  which  is  in  sharp  contrast  with  the 
assumption  of  an  isotropic  distribution  of  electrons  in  the 
momentum  space  adopted  by  the  energy-balance  equation.3 
Relatively  cool  electrons  in  kx> 0  states  (sharp  tail  at  the 
Fermi  surface)  are  found  in  comparison  with  hot  electrons  in 
kx< 0  states  (smooth  tail  at  the  Fermi  surface).  Moreover,  we 
see  by  comparing  Fig.  8(b)  with  Fig.  8(a)  that  the  depen¬ 
dence  of  f(sk  | ,  ek )  on  k  and  kx  is  also  anisotropic,  and  the 
electrons  are  relatively  unheated  or  cooled  (almost  identical 
Fermi  surface)  by  the  dc  field  in  the  direction  perpendicular 
to  the  field.  This  excludes  the  possibility  of  defining  a  state- 
independent  electron  temperature  as  introduced  by  the 
energy-balance  equation.3 

Figure  9  illustrates  the  field  dependence  of  calculated  drift 
velocity  u0  in  Fig.  9(a)  and  mobility  /j.c  =  uq/ Edc  in  Fig.  9(b) 
at  T=30  K  when  the  energy-drift  term  in  Eq.  (10)  is  included 
(solid  curves  with  the  symbol  □)  and  excluded  (dashed- 
dotted  curves  with  the  symbol  O).  In  Fig.  9(a)  we  can  see 
that  the  energy-drift  effect  heats  the  electrons  in  the  field 
direction  and  reduces  the  drift  velocity  when  £dc  is  strong 
(greater  than  0.75  kV/cm).  The  nonlinear  relation  between 
the  drift  velocity  and  the  applied  dc  field  (solid  curve)  can  be 
seen  very  clearly  in  the  strong-field  regime  when  the  energy- 
drift  effect  is  included,  which  is  an  indication  of  the  high- 
field  transport  of  hot  electrons.  Moreover,  we  find  in  Fig. 
9(b)  that  the  mobility  /xt  decreases  with  EAc,  which  has  been 
observed  previously  from  numerical  calculations,10  and  the 
energy-drift  effect  further  reduces  the  mobility  when  Edc  is 
higher  than  0.75  kV/cm. 

Finally,  we  conclude  from  Figs.  8  and  9  that  the  energy- 
drift  equation  derived  in  this  work  does  include  the  addi¬ 
tional  hot-electron  effect  in  the  high-field  transport  in  com¬ 
parison  with  the  regular  Boltzmann  transport  equation. 
Furthermore,  the  energy-drift  equation  also  includes  the  an¬ 
isotropic  distribution  of  electrons  in  momentum  space  along 
and  perpendicular  to  the  field  direction,  due  to  the  accelera¬ 
tion  of  electrons  by  the  dc  electric  field.  This  implies  that  the 


FIG.  8.  Comparison  of  calculated  dimensionless  electron- 
distribution  functions  f(Ekx,ek)  =  (sf/3TT1(Tw)f(ek^,sk)  as  func¬ 
tions  of  Ek  =Ek  / Ef  and  Ek^  =  s kJ £f  in  (a)  and  (b)  at  T=  30  K  and 
£dc  =  0  .75  kV/cm  for  hot  and  equilibrium  electrons.  The  dashed- 
dotted  curves  in  both  (a)  and  (b)  with  the  symbol  □  represent  the 
results  from  the  equilibrium  electrons.  In  (a),  k±  =  0  and  the  solid 
and  dashed  curves  with  symbols  O  and  A  represent  the  results  from 
the  hot  electrons  with  kx< 0  and  kx> 0,  respectively.  In  (b),  kx= 0, 
and  the  solid  curve  with  the  symbol  O  represents  the  result  from  the 
hot  electrons. 

isotropic  assumption  adopted  by  the  energy-balance  equation 
cannot  be  justified  under  a  strong  dc  electric  field.  On  the 
other  hand,  when  the  dc  field  is  weak,  our  coupled  energy- 
drift  and  force-balance  equations  reduce  to  the  same  results 
as  those  obtained  from  coupled  energy-balance  and  force- 
balance  equations,  as  well  as  from  the  regular  Boltzmann 
transport  equation,  when  only  the  linear-field  terms  are 
kept.3,4 


IV.  CONCLUSIONS  AND  REMARKS 

In  conclusion,  coupled  energy-drift  and  force-balance 
equations  have  been  derived  for  high-field  hot-electron  trans- 
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ation  in  light-emitting  diodes  or  quantum-well  lasers,11  the 
study  of  saturation  drift  velocities  in  different 
semiconductors,12  the  study  of  phonon-drag  thermoelectric 
power,9  and  the  study  of  the  nonresonant  interaction  of  car¬ 
riers  with  terahertz  radiation. 13 

By  further  incorporating  the  electric  and  magnetic  poten¬ 
tials  into  the  proposed  energy-drift  and  force-balance  equa¬ 
tions,  we  will  be  able  to  explain  nonlinear  miniband  trans¬ 
port  in  superlattices14  and  the  magneto-transport  effect  in  a 
two-dimensional  electron  gas  (2DEG).15  Moreover,  by  in¬ 
cluding  the  energy-drift  term  in  the  semiconductor  Bloch 
equations,7  we  will  be  able  to  explain  transport  effects  on 
laser-induced  optical  coherence  and  vice  versa. 
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APPENDIX 


By  using  the  three  dimensionless  functions  defined  in 
Eqs.  (37)-(39),  we  can  express  all  the  dimensionless  expan¬ 
sion  coefficients  of  the  Fokker-Planck-type  equation  as 


^T(e*±>  ek) 


hAT(ski,sk) 

Ef 


1  f+0° 

— S  dqLq± 

lit  X  JO 


FIG.  9.  Comparison  of  calculated  drift  velocities  u0  and  mobili¬ 
ties  fic  =  u0/Edc  as  functions  of  the  electric  field  Eic  in  (a)  and  (b)  at 
T=30  K  with  (solid  curves  with  the  symbol  □)  and  without  (dashed 
curves  with  the  symbol  O)  an  energy-drift  term  in  the  energy-drift 
equation. 
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r+ co 
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J  —00 


1  r+ oo  r+ 00 

—  dq±q±  dqx\Cq\\itr+{{k},{q}) 
2'77-Jo  J_m 


port.  The  work  done  by  the  frictional  force  has  been  included 
in  the  Boltzmann  scattering  equation  for  electron  relative¬ 
scattering  motion  and  has  been  found  to  increase  the  thermal 
energy  of  the  electrons  and  to  reduce  the  field-dependent 
drift  velocity  and  mobility  of  the  electrons.  The  importance 
of  the  hot-electron  effect  in  the  energy-drift  term  under  a 
strong  dc  field  has  been  demonstrated.  The  Doppler  shift  in 
the  energy  conservation  of  scattering  electrons  interacting 
with  impurities  and  phonons  has  been  found  to  give  rise  to 
the  anisotropic  electron  distribution  in  momentum  space 
along  the  field  direction.  The  importance  of  this  anisotropic 
distribution  has  been  demonstrated  through  a  comparison 
with  the  isotropic  energy-balance  equation,  from  which  the 
possibility  for  defining  a  state-independent  electron  tempera¬ 
ture  has  been  excluded. 

The  proposed  energy-drift  and  force-balance  equations 
can  provide  new  physical  insight  in  several  areas.  These  in¬ 
clude  the  study  of  current-injected  hot-carrier  energy  relax- 
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(A5) 


where  n,  =  Na/V  is  the  concentration  of  impurity  atoms,  sk 
=  skJsf,  Bk=skJsf,  Ui(q)  =  (Ze2kf/e0ersf)/(q2±  +  q2x  +  Q2), 

Q2  =  (e2kjl e0erSf)l (2tt).  Moreover,  we  find  the  dimension¬ 
less  form  for  the  electron-phonon  coupling. 


|CJ2  = 


_  e2kfill  1  j_ 


£o  ef 


1 


Cm  £<■  /  tP-  j.  77-  j.  zi2 


(A6) 


<7i  +  C 


|C«P=l^ 


2piSf/  . 


D2(q2±  +  qx)  +  32(ehu)2 


(q\  +  q2xf2 


J  (ql  +  ql+Q2)2 
(Al) 


\2pis,/64  (o2  +  o?  +  ( 


,2x2 


(<7l  +  lx  +  2j) 


(A8) 


where  D=D/ Sp  ehl4=ehl4/ kj-ep 

By  using  the  dimensionless  functions  defined  in  Eqs. 
(37)-(39),  we  can  express  all  the  dimensionless  emission 
and  absorption  rates  of  the  dynamical  equation  for  hot 
phonons  as 


h@T 

em _ 


@-m  = 


r+ 00  r 

=  2tt|C^|  2J  dkLkL^ 


dk±k±  I  dkx(k±  +  k~ ) 


Z 2  ,  E2X-1/2 


x  [&({k},{q})  +  €\({^K{<7})]  +  ^'ACifh\1{qsk  -  qxu0) 


r+oo  r  4 

s  I  dk±k±  I 

Jo  J  - 


x  I  dkLkL  I  ^,(P±  +  ^)-1/2^(W,®) 


/  niU20  \ 

P  P 

\  3  7^0-30  / 

X 

1 

H 

+ 

^  1 

/'I 

1 _ 

/(e*:±,e*r), 


(A9) 
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r+x>  r+ 00 

•/  0  — 00 

2«oV  r\_  _ 

- 2j  dqLq± 

tt  x  Jo 

r+  00 

dqxqx(qxu0-qs\)\Cqx\2 

J  —CO 


abs _  #A 


/j^aos  _ 

U<A  - 


=  2tt|C?-x|2 


r+co  r 

J  PM,  J 


dk±k4  I  dkx(kl±  +  kx) 


E2X-1/2 


X  [^(W,{<?})  +  £I(W,  {<?})]  -  47r|C^|2(^x  -  qxu0) 


X{[M-x(2^x)  +  l]®},®) 
+  ^(2^)  £({£},  ®)} 


r+co  r  + 

'  I  dkJiL  J 
Jo  J-c 


x|  ^i+^)-i/2^(w,{?}) 


X 


p  p 

Z  Z 

dek  dek 


f(Ek±’ek)- 


(A10) 


2«0  r+o°  _  _  r+c _  - .  - . . 

Jo  J  -00 


For  optical  phonons,  we  have  similar  results. 
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Sb 
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+ 
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(All) 

X 
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1 1  is 

1  ^ 

+ 
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f(skx,sk)- 
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